Nonlinearity in Bacterial Population Dynamics: 
Proposal for Experiments for the Observation of Abrupt Transitions in Patches 
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An explicit proposal for experiments leading to abrupt transitions in spatially extended bacterial 
populations in a Petri dish is presented on the basis of an exact formula obtained through an analytic 
theory. The theory provides accurately the transition expressions in spite of the fact that the actual 
solutions, which involve strong nonlinearity, are inaccessible to it. The analytic expressions are 
verified through numerical solutions of the relevant nonlinear equation. The experimental set-up 
suggested uses opaque masks in a Petri dish bathed in ultraviolet radiation as in Lin et al., Biophys. 
J. 87, 75 (2004) and Perry, J. R. Soc. Interface 2, 379 (2005) but is based on the interplay of two 
distances the bacteria must traverse, one of them favorable and the other adverse. As a result of this 
interplay feature, the experiments proposed introduce highly enhanced reliability in interpretation 
of observations and in the potential for extraction of system parameters. 



Phenomena displaying abrupt transitions are of special 
interest to a variety of sciences, including physics and bi- 
ology. In some physical cases they arise from cooperative 
interactions among a large number of constituents (e.g., 
molecules or spins), and in others from nonlinearities in 
interaction inherent in the system, (e.g., in some mechan- 
ical systems with only a few degrees of freedom^ In the 
first form they are known as phase transitions [lj , in the 
second as bifurcations [2j]. Abrupt phenomena also com- 
mand attention in the context of extinction of popula- 
tions, a subject of obvious interest to biology Q. Tran- 
sitions in populations, therefore, constitute an exciting 
topic of interdisciplinary science combining physics and 
biology, and the present paper reports a theory and pro- 
poses an experiment in this topic. The special feature 
of the proposed experiment is that it may be performed 
with relatively simple equipment and measurement tech- 
niques. 

The system to be considered consists of bacteria in 
a Petri dish, allowed to grow and move, spatial selec- 
tion being imposed via lethal ultraviolet radiation that 
is incident on the dish but punctuated by opaque masks 
that protect the bacteria in chosen regions. Experiments 
under such a setup were initiated several years ago by 
Lin et al. [3| who used moving masks in response to a 
theoretical analysis [5j that focused on spatial disorder. 
A quite different experiment using stationary masks of 
varying sizes, and employing bacterial extinction as the 
key phenomenon was then proposed @. That proposal 
was carried out experimentally 0] and the predictions of 
ref. [f| were verified. Related investigations on this topic 
may be found in [j| 0, |T(| , the general problem of bac- 
terial dynamics in Petri dishes having been addressed in 
various papers and contexts 

The underlying assumption behind most of these stud- 
ies is that the bacteria obey a simple Fisher equation [l2j 
for the time evolution of their dynamics, their population 
density u(x, t), where x is the position in a 1-dimcnsional 
space (the linear dimension of the Petri dish) , and t is the 
time, being governed by 



du{x) 
dt 



= D 



d 2 u{x) 
dx 2 



+ au{x) — bu 2 (x). 



(1) 



Here a and b are, respectively, the growth rate and a 
competition parameter arising from the resources being 
limited, and D is the bacterial diffusion constant. 

Equation |T]) generally does not permit analytic solu- 
tions for arbitrary times but can be solved in the steady 
state explicitly in terms of Jacobian elliptic functions, 
if the bacterial population is assumed to vanish at the 
edges of a finite region. Physically, this could correspond 
to an idealization in which, while the mask protects the 
bacteria from the lethal effects of the ultraviolet radi- 
ation, the latter is so potent that bacterial population 
must vanish everywhere outside the mask where the ra- 
diation impinges on the bacteria. Mathematically, these 
are Dirichlet boundary conditions. The elliptic function 
solution, given independently in a recent analysis [|| , but 
known decades earlier through the work of Skellam [l3j |. 
and elucidated in textbooks [15j , leads to the well-known 
KISS transition [lij ]: the steady state bacterial popula- 
tion vanishes for any width of the mask lower than a 
critical value that, interestingly, depends on the diffusion 
constant D and the growth rate a, but not on the nonlin- 
earity parameter b. The specific expression for the crit- 
ical mask width is %\J D/a. Combining this expression 
with information about bacterial diffusion constants and 
growth rates discussed by Mann in his Ph.D. thesis 
Kenkre and Kuperman calculated the critical mask width 
to be around 0.5cm, and suggested that an experiment be 
carried out to observe the transition. Perry Q followed 
the suggestion, and reported observing a transition width 
of 0.8 cm in his experiments on a non-chemotactic strain 
RP9535 of E. coli bacteria. Since the Dirichlet boundary 
condition is only an extreme idealization, Perry argued, 
appropriately, that it is preferable to invoke the anal- 
ysis of Ludwig et al. [161 ] for quantitative verification. 
That analysis assumes that the growth rate a is nega- 
tive in the region outside the mask (and of magnitude 
a\) rather than infinite as would correspond to Dirich- 
let boundary conditions, and arrives at the critical mask 
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FIG. 1: Our proposed experimental set-up showing two fa- 
vorable regions each of length L separated by a hostile region 
of length 2d. The focus of observation is the interdependence 
of the critical values of L and d at which extinction occurs. 

width as being L = 2 \J~Dja arctan a\j a. The Dirich- 
let expression is recovered from this Ludwig formula if 
the effect of the ultraviolet radiation outside the mask is 
"infinitely" lethal. 



I. PROPOSAL FOR EXPERIMENT 

Need for the new experimental proposal that we 
present here arises from the fact that the stationary sin- 
gle mask experiment @, @] uses for its interpretation a 
number of inputs which suffer from a certain degree of 
uncertainty. Some of these uncertainties are about the 
values of D and a [9(, and others about the extent to 
which other phenomena such as signaling [7] that occur 
in bacterial movement might affect the outcome of the 
experiment. Because absolute values of critical lengths 
may be difficult to obtain with an acceptable degree of 
accuracy, we propose here stationary mask observations, 
that focus, instead, on the interdependence of two crit- 
ical lengths. The essential feature of this present pro- 
posal is to have more than one mask in the Petri dish, 
so that there are two controllable lengths, a favorable 
length associated with a livable region for the bacteria, 
and a hostile length associated with an unlivable region. 
The double control possible in the proposed experiment 
opens the prediction space from a single value to an in- 
finity of values. The variation of the critical value of one 
of the lengths when the other is changed presents an en- 
tire relationship that is directly observable rather than a 
single value and therefore could lead to a much cleaner 
and trustworthy interpretation of the observations. 

The simplest system to study consists of two masks, 
each of length L, separated by a distance of length 2d as 
shown in Fig. [H The previously investigated case with 
a single mask [(| @] corresponds, obviously, to infinite d. 
To test the idea behind the new proposal, we carried out 
numerical studies of the full nonlinear Fisher equation 
([T|) with given values of D, a, a\, and b (10, 0.1, -0.9 
and 1 in appropriate units.) All lengths were expressed 
as ratios to the diffusion length y^D/a. The numerical 
studies used an explicit finite differences scheme, x be- 
ing discretized in intervals of 0.1. The convergence to 
a steady state was analyzed by measuring the distance 
between successive solutions. 

The result we found ( see Fig. [2] ) is that the critical 
value of the mask width L, i.e., the smallest value that 
can support a non-zero bacterial population in the steady 




FIG. 2: Basis of the experimental proposal: Shown is the 
interdependence of the critical values of the favorable and 
hostile lengths L and d respectively. The former is the width 
of the mask and the 2d is the distance between the two masks 
(see Fig. 1), both lengths being expressed in this plot in 
units of the diffusion length ^ D/a where a is the growth 
rate under the masks. Outside the masks a is replaced by 
— ai to represent the harsh effect of the ultraviolet light, the 
magnitude of a\ being taken for the purposes of this plot to 
be 9 times that of a. We solved Eq. [T] numerically for various 
different values of the nonlinear parameter b and found that 
the results did not depend on those values provided b was non- 
zero and positive. The shaded region represents pairs of L and 
d values that lead to bacterial extinction in the steady state. 
Filled circles mark the onset of extinction and are obtained 
from the numerical solution. The solid line constructed simply 
to smoothly join the circles, is found to coincide precisely with 
the prediction of our analytic theory. 



state, is lowered for a finite intermask distance d. We con- 
sidered 10 different values of the hostile length (the in- 
termask distance d) and varied the favorable length (the 
mask width L.) For each L-d pair, we started with arbi- 
trary initial conditions and let the program run until no 
time dependence was discernible. 

We repeated the procedure for each of several suffi- 
ciently low values of L and increased L until the extinc- 
tion disappeared. We also reversed the procedure start- 
ing with high values of L and decreased them systemat- 
ically until extinction appeared. Numerous runs allowed 
us to obtain corresponding pairs of L and d that mark the 
transition region. The results are denoted by filled cir- 
cles in Fig. [2] The shaded area represents the extinction 
region and the unshaded area the parameter region in 
which bacterial population densities are non-zero in the 
steady state. The curve passing through the numerically 
found transition points may be considered, at this stage 
of our discussion, to be simply a smooth joining trace. 
We will see below that its exact shape can be accessed 
through our analytic theory. 

It is easy to understand, on the basis of a qualitative 
argument, the shape and tendency of the results of the 
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numerical solutions of the nonlinear equation. Bacteria 
diffuse from within the mask to the harsh region and 
die if they reach that region. Small values of the mask 
width or large values of the intermask distance result 
in extinction. The extinction effect is worsened by an 
increase in the intermask distance for small values of that 
distance but the effect saturates for larger values. Hence 
the saturation in the curve. 

All these features can be tested experimentally in our 
proposed set-up. Quantitative comparison with the pre- 
dictions of the Fisher equation are possible because we 
have developed an exact analytic theory of the inter- 
play of the favorable and hostile distances as reflected 
in the transition. What makes the proposed comparison 
with experiment significant is that, although the Fisher 
equation, whose numerical solution has led us to set out 
the separation curve between the extinction region and 
the rest in Fig. [T] cannot be solved exactly by analytic 
means, the separation curve itself can be obtained ana- 
lytically. Indeed, we show below that the curve is given 
precisely by 




FIG. 3: Our analytic predictions for the dependence of crit- 
ical values of the favorable length L and the hostile length 
d respectively from Eq. @. As in other plots, the lengths 
are expressed as dimensionless ratios to the diffusion length 
sj D la. Different curves correspond, as shown in the legend, 
to different values of £, i.e., to different intensities of the ul- 
traviolet radiation. Note the doubling of L as d passes from 
to oo in all curves. 



L = \ — 



arctan w h arctan 



— tanh 




(2) 



This prediction, which is one of the central analytic re- 
sults of the present paper, coincides with the solid curve 
in Fig. EJ 



inside the mask, d < x < d + L, 

u(x) = O exp(— a.\ \x\) 



(0) 



in the extreme outside, in the harsh region, x > d + L, 
and 



II. ANALYTIC THEORY FOR THE TWIN 
MASK SET-UP 

Following the ideas of Ludwig et al. [l6[ but applying 
them to the many-mask system, we consider the steady 
state of the Fisher equation |T|) and argue that if there is 
a transition, the quadratic term in the steady state u{x) 
can be neglected at the extinction point in favor of the 
linear terms since u(x) itself vanishes at the transition. 
We are thus led to seek the solutions of 



d 2 u(x) 
dx 2 



+ a 2 u(x) = 



under the two masks, and of 



d 2 u(x) 
dx 2 



a\u(x) = 



(3) 



(4) 



outside the masks, with a 2 — a/D and a\ = a\/D. We 
take the masks to lie from x = ±d to x = ±(d + L). 
Using / and O as constants, the most general functions 
for u{x) inside and outside the mask, respectively, are, as 
a result of the symmetry, (we consider only the right side 
of the origin, since all considerations repeat unchanged 
on the left side by symmetry) 



u(x) — /cos(a|x| 



(5) 



u(x) — Ccosh(ai|a;|) 



(7) 



in the central region between the two masks, —d < x < d. 
Matching the logarithmic derivative of the solution at the 
outer and the inner boundaries leads to 



tan[a(d + L) — 
tan(arf — < 



= cti/a 

= — (ai/a) tanh ct\d, 



and elimination of <j) from these two equations leads to 
Eq. ([2]) quoted above. Notice how, in light of the behav- 
ior of the hyperbolic tangent, Eq. @, one of our central 
results, reduces to Ludwig et al.'s single mask value [l6| 
of L, when d attains infinite values, and half that value 
when d vanishes. Both are fully expected and natural 
results. It is also instructive to rewrite the transition 
relation as 



C = arctan £ + arctan(£ tanh 5£) 



(8) 



where we express the favorable distance L and hostile 
distance d normalized to the growth diffusion length g = 
\f~D~ja~ as C — L/g and 5 — d/g respectively, and the 
depletion parameter £ is the square root of the ratio of 
the destruction rate a\ and the growth rate a. 

We display in Fig. [3] our analytic result, Eq. |(2]), 
equivalently Eq. ([8|). As in other plots the lengths are 
expressed as dimensionless ratios to the diffusion length 

D /a. Different curves correspond, as shown in the leg- 
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end, to different values of £, i.e., to different intensities 
of the ultraviolet radiation. The curve for £ = 3 from 
our theory is shown as the solid line in Fig. [21 and dis- 
plays exact coincidence with the numerical findings. To 
be noted is the fact that the numerical solutions are of 
the full nonlinear equation and have ben found for a spe- 
cific b whereas the analytic theory is linear and does not 
require the use of a value of b, requiring for its application 
only the condition b > 0. 

It is clear from our Eq. ((2| that the variation of the 
destructive rate a\ may also be employed for useful ex- 
perimental exploration, at least in principle. While this 
could have also been done in the one-mask scenario of 
ref. @, we explain the idea here in the two-mask case. 
Additional confirmation of the theoretical picture of the 
dynamics of the bacteria might be provided by observing 
how the critical L-d values change with the intensity of 
the ultraviolet light. Even if we do not know the precise 
dependence of a\ on the intensity, we can be fairly certain 
that it increases with the latter, and undergoes a satu- 
ration at high values of the intensity. A variation of the 
critical value of L with the intensity would show there- 
fore qualitative behavior similar to the variation with a\ , 
equivalently with the ratio £ = y a\/a. The latter vari- 
ation is displayed in Fig. 0J With one exception, all d 
values produce the same saturation value of L for large 
£. The exception is d = 0: its saturation value is one half 
that of the others. Mathematically, this corresponds to 
the hyperbolic tangent becoming 1 for all non-zero val- 
ues of d for large enough £, but vanishing if d vanishes. 
Physically, this means that if there is an adverse region 
between the masks, bacteria will be killed on arriving 
there, in light of the infinitely harsh radiation, reducing 
the problem to a single-mask scenario with the given L as 
the mask width. On the other hand, if the intermediate 
harsh region does not exist at all (because d — 0), one 
is reduced to considering a single mask with width twice 
that of the given mask. 

Our suggestion, thus, is to use both the relative varia- 
tion of the favorable and hostile lengths, L and d, and of 
the intensity of ultraviolet light in the manner discussed, 
to check our simple quantitative predictions. 



III. MULTIPLE MASK SET-UP 
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FIG. 4: Our analytical prediction from Eq. for the depen- 
dence of the critical value of L on the intensity of ultraviolet 
radiation, represented here by the value of £, for given val- 
ues of the hostile length (half the intermask separation) d as 
shown. All d values produce the same saturation values of 
L for infinitely harsh ultraviolet intensity except for d — 0: 
all curves reach the same higher value for large £ except the 
d = curve whose saturation value is one half that of the 
others. See text for explanation. The lengths a re ex pressed 
as dimensionless ratios to the diffusion length yjD/a. 



angent of the result, and (iv) subtracting the result from 
£. Define the operator T similarly through the successive 
operations of (i) taking the trigonometric tangent of what 
it acts on, (ii) dividing the result by £, (iii) taking the 
hyperbolic arctangent of the result, and (iv) subtracting 
the result from 2£<5.To obtain the required formula for an 
even number of masks, start with the product £<5. Apply 
H and then T alternately and successively so that there 
are as many £Ts in the operation as the number of mask 
pairs in the system, but one more H than T. (Thus, do 
H for a pair of masks, HTH for 4 masks, HTHTH for 
6 masks, and so on.) Finally, equate the result to the 
arctangent of £ to get an implicit formula showing the C- 
5 relation. The pseudocode can be expressed succinctly 
in terms of operators T and H by stating that the gen- 
eral expression for 2n masks can be formally written as 



(£5) = arctan £. As can be verified, this 



reproduces the 2-mask formula, 
formula is 



Eq. ©. The 4- mask 



It is possible and useful to construct set-ups with mul- 
tiple masks for further experimental verification. The 
theory for this situation along the Ludwig et al. argu- 
ments we have developed here is slightly more tedious to 
write down. For ease in notation, let us adopt as we have 
done in Eq. (jBJ), the symbols C and 5 to express the fa- 
vorable and hostile lengths in units of y/D/a. Then the 
linear Petri dish formula for an even number of masks is 
obtained by executing the following pseudocode. 

Let H be an operator defined by the successive opera- 
tions of (i) taking the hyperbolic tangent of what it acts 
on, (ii) multiplying the result by £, (iii) taking the arct- 



C = arctan £ + arctan ( £ tanh (2£<5-arctanh((I/£) 
tan(£-arctan(£tanh £<5))))^. (9) 

IV. CIRCULAR PETRI DISH 

Our suggested set-up and calculations have taken the 
Petri dish to have infinite extension away from the two 
masks, an assumption that should be reasonable in light 
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FIG. 5: The proposed set-up in a circular geometry shown 
here with a multiple mask set-up with 3 experimentally con- 
trollable parameters: R as well as L and d. 



of the destructive effects of the radiation. We now ana- 
lyze another practical extension of the theory. Consider 
the two-mask situation in a circular Petri dish, i.e. one 
in which the boundary conditions are periodic, the hos- 
tile distance between the edges of the masks being 2d on 
one side as in the earlier analysis but now also 2R on the 
other side, see Fig. [5] Let the width of the dish be small 
enough so that bacterial diffusion can be considered to 
be still 1-dimensional but over a total distance of extent 
2(L + d + R). We are not interested here in 2-dimensional 
considerations appropriate to wide dishes, an example of 
which may be found in ref. [lj} for experiments with 
moving masks or rotating dishes. 

We now have, as in the analysis given above, u(x) = 
Jcos(a|x| — 4>) inside the right mask, i.e., for d < x < 
d+L, but to the right of the right mask, i.e., for x > d+L, 
u(x) is not proportional to exp(— ai|x|) as in the earlier 
analysis, but 



u(x) = O cosh(ai \x\ — rf). 



(10) 



This expression, which, in the right region under con- 
sideration, where x > 0, can be written as u(x) — 
O cosh(aia; — 77), will hold on the left extreme as well but 
there it simplifies to u{x) = Ocosh(aia; + rf). The addi- 
tional phase 77 is immediately determined from the peri- 
odicity in the boundary conditions, i.e., by matching the 
logarithmic derivative of the solution at \x\ = d + L + R: 



i] = a>i(d + L + R). 



(11) 



Petri dish: 



L = 




(12) 



We see that the circular dish result, Eq. (fT2")l , is sym- 
metric in the two hostile lengths R and d as it obviously 
must be, that it reduces to the result for the linear dish 
of infinite extent, viz., Eq. when R becomes infinite 
since then the new hyperbolic tangent in Eq. (112p reduces 
to 1, and that, when R vanishes (equivalently when d van- 
ishes), it yields a new expression for the critical width for 
the single-mask case but with twice the width value. If 
there were simply a single mask in a circular Petri dish 
of extent R and no d, the generalization of the Ludwig 
formula that would come out of our analysis would be 




arctan 



tanh(aii?) 



(13) 



Most importantly, we see that it should be possible in 
principle to use the new analytic result (s) for experimen- 
tal probing of the phenomenon under consideration since 
one could construct and employ circular Petri dishes of 
desirable radius. Needless to say, the added effects dis- 
cussed here would not be discernible for R values that 
are very large as then bacteria would die in the external 
regions quickly enough. 

We have also derived the corresponding expressions 
combining circular dish geometry with the multiple mask 
set-ups. Specifically, the 4-mask fo rmul a takes the form, 
1Z being R measured in units of y/D/a, 



C = arctan (£tanh£7?.) + arctan f^tanh ^2£<5 — arctanh( 
(l/£)tan(£ - arctan(£tanh ££))))). (14) 



Matching also at x = ±d and x — ±(d+L), one arrives at 
the new formula for critical L and d, valid for a circular 



whose reduction to the linear dish formula is obvious as 
R becomes infinite. We have verified all these expressions 
that we have derived, by comparing their predictions to 
numerical solutions of the Fisher equation with explicit 
(arbitrary) non-zero 6's. 



V. CONCLUDING REMARKS 

The proposal for experimental observations of abrupt 
population transitions in patches that we have presented 
above should be of interest for multiple reasons. Patches 
or spatial inhomogeneities have been studied and empha- 
sized early on in ecology by many authors such as Levin 
[III, Murray [l9| and Shigesada [HJ in varied contexts 
including traveling waves, and there have been recent 
contributions as well on the Dirichlet conditions problem 
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for general nonlinearity incorporating Allee effects [211 ]. 

Our analysis here provides a clean proposal for ex- 
periments which derives its potential for clear interpre- 
tation from the use of two (or more) controllable dis- 
tances of travel, one favorable and the other unfavorable. 
One should be able to use it for extraction of parame- 
ter combinations in the Fisher equation which has been 
ubiquitous in mathematical ecology. The idea behind 
the present proposal should therefore find use in systems 
other than bacterial aggregates. Indeed, we have found 
similar abrupt transitions in infected rodent populations 
in our study of the Hantavirus epidemic. The theory 
for that system is, however, more complex, and cannot 
be developed without approximation. By contrast, the 
theory we have presented here is, perhaps surprisingly, 
exact, i.e., analytical. This is so in spite of the fact that 
analytic solutions of the Fisher equation for these situa- 
tions do not seem to be possible except as elliptic quadra- 
tures. The reason for this happy state of affairs can be 
understood from the original arguments of Ludwig et al. 
given in their paper on spatial patterning of the bud- 
worm [16], or from the lucid explanations given, e.g., in 
a recent text on mathematical ecology [HI . At the tran- 
sition, the densities vanish and therefore terms of order 
higher than the first may be safely neglected. Our theo- 
retical contribution is only in generalizing that analysis 
to multiple masks and nonlinear geometries as shown. 
The observationable controllable features our proposal 
emphasizes, and directly utilizes, are (i) the interplay of 



the favorable and adverse lengths for bacterial traversal, 

(ii) the variation of the intensity of ultraviolet radiation 
which changes the interrelationship of the critical values, 

(iii) the multiple-mask set-up which reintroduces bacteria 
into favorable regions after they have passed into arid ar- 
eas, and (iv) the optional use of circular geometry for the 
Petri dish which provides one more controllable length. 

We are currently involved in investigating fluctuation 
effects that appear from detailed Monte Carlo consider- 
ations and become manifest for small population densi- 
ties, in constructing the theory for transitions in partially 
infected populations of rodents in open terrains, in ana- 
lyzing the effects of static disorder in the placement and 
size of the masks, and in studying the effects of dynamic 
(time) variation of mask size and placement as in ref [Io| . 
The bacterial experiments proposed above are, however, 
ready to go and we hope that observations will be made 
soon along the lines we have discussed. 
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